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CN ; Abstract 

O . We justify supercritical geometric optics in small time for the defocusing semiclassical Non- 

.^ ' linear Schrodinger Equation for a large class of non-necessarily homogeneous nonlinearities. The 

case of a half-space with Neumann boundary condition is also studied. 
OO . 

PLi ■ 1 Introduction 

<'■ 

• ' We consider the nonlinear Schrodinger equation in fi C M^ 

a' ie^ + — A*^-^V(|^^|2) = 0, ^':R+xn^C (1) 

with an highly oscillating initial datum under the form 

in: ^^, = % = alexp(^^^'^, (2) 

(N' 

T-H ' where c/9g is real-valued. We are interested in the semiclassical limit e ^ 0. The nonlinear 

Ti^lj- . Schrodinger equation ([1]) appears, for instance, in optics, and also as a model for Bose-Einstein 

^D I condensates, with f{p) = p — 1, and the equation is termed Gross-Pitaevskii equation, or also 

with f{p) = p^ (see |13j). Some more complicated nonlinearities are also used especially in low 

dimensions, see [12] . 

.^H ! At first, let us focus on the case i7 = M.'^. To guess the formal limit, when e goes to zero, it is 

kS ' classical to use the Madelung transform, i.e to seek for a solution of ([T]) under the form 

^^ = ^/Jfexp ( -93^ 

By separating real and imaginary parts an by introducing u^ = Vc/5^, this allows to rewrite ([T]) as 
an hydrodynamical system 

' dtp"" + V • {p^u') = 

dtu^ + {u' ■ vy + v(/(p^)) = -V ( -^ ) • 
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The system ([3]) is a compressible Euler equation with an additional term in the right-hand side 
called quantum pressure. As e tends to 0, the quantum pressure is formally negligible and ([3]) 
reduces to the (compressible) Euler equation 



dtp + V- {pu) = 
dtu+{u-V)u + V{f{p))=0. 



(4) 



The justification of this formal computation has received much interest recently. The case of analytic 
data was solved in [7] . Then for data with Sobolev regularity and a defocusing nonlinearity, so that 
@ is hyperbolic, it was noticed by Grenier, [9], that it is more convenient to use the transformation 



^' 



exp [i — 1 



(5) 



and to allow the amplitude a^ to be complex. By using an identification between C and M^, this 
allows to rewrite (fl]) as 



f a"- 

dta^ + u"" • Va^ + — V-u^ = eJ Aa^ 



(6) 



^ dtu^ + {u^ ■ V)u^ + V(/(|a=p)) = 0, 
where J is the matrix of complex multiplication by i: 



J 



-1 

1 



When e = 0, we find the system 



dta + u ■ Va + - V • it = 
_ dtu+{u-V)u + V{f{\a\^)) =0, 



(7) 



which is another form of dU, since then (p = japju) solves ([!]). The rigorous convergence of ^ 
towards ([7|) provided the initial conditions suitably converge was rigorously performed by Grenier 
[S] in the case f{p) = p (which corresponds to the cubic defocusing NLS). More precisely, it was 
proven in [9] that there exists T > independent of e such that the solution of ([6]) is uniformly 
bounded in H^ on [0,T]. In terms of the unknown ^^ of ([1]), this gives that 



sup sup II ^^ exp 

ee(0,l] [0,T] 



■'P^ 



\H' 



< +00 



for every s where {a,u = Vc^) is the solution of ([7]). Furthermore, the justification of WKB 
expansions under the form 



^' 



e a 



e e 



Oie' 



k=0 



for every m was performed in [9]. The main idea in the work of Grenier [9] is to use the symmetrizer 

S^ diag(l,l,i/'(|a|2),...,i/'(|a|2; 



of the hyperbolic system ([7]) to get H^ energy estimates which are miiform in e for the singularly 
perturbed system Q. The case of nonlinearities for which /' vanishes at zero (for instance the 
case /(p) = p"^) was left opened in |9j. The additional difficulty is that for such nonlinearities, the 
system ([7]) is only weakly hyperbolic at a = and in particular the symmetrizer S is not anymore 
positive definite at a = 0. 

In more recent works, see |19j . |14j . [l] it was proven that for every weak solution of ([T]) with 
f{p) = p — 1 or f{p) = p, the limits as e — > 

|\&=|2-p^0 in L°°([0,T],l2) elm {^^V-^^) - pu ^ in L'^{[0,T],lIJ (8) 

hold under some suitable assumption on the initial data. The approach used in these papers is 
completely different, and relies on the modulated energy method introduced in [3]. The advantage 
of this powerfull approach is that it allows to describe the limit of weak solutions and to handle 
general nonlinearities once the existence of a global weak solution in the energy space for ([T]) 
is known. Nevertheless, it does not give precise qualitative information on the solution of ([1]), 
for example, it does not allow to prove that the solution remains smooth on an interval of time 
independent of e if the initial data are smooth or to justify WKB expansion up to arbitrary orders 
in smooth norms. 

In the work [2], the possibility of getting the same result as in [9j for pure power nonlinearities 
f{p) = p" in the case ri = M'^ was studied. It was first noticed that, thanks to the result of |15j . 
the system 

dta + V(/3 • Va + -Ay? = 

1 (9) 

d,^ + -\V^\^ + f{\a\^) =0, 

with the initial condition {a,(f^ , = (ao,(/3o) £ H'^ has a unique smooth maximal solution 
(a,v3) G C([0,r*[,i7^(M'^) X H^-'^{W^)) for every s. It was then established: 

Theorem 1 (^) Let d < 3, cr G N* and initial data a^, c/Jq i'n- H°° such that, for some functions 

ho-(^o\\Hs = C{e), 

for every s > 0. Then, there exists T* > such that ([9]) with f{p) = /9°" has a smooth maximal 
solution {a,(p) € C([0,T*[, i7°° x H°°). Moreover, there exists T £ (0, T*) independent of e, such 
that the solution of ([T|), d^J remains smooth on [0,T] and verifies the estimate 



' ^J\\l°°{[o,t],h^) 



sup ||*^exp(-z-)||^<^^^Qy^^,-| < +00, (10) 

£6(0,1] 



where 

• if a = 1, then s £ N is arbitrary, 

• if a = 2 and d = 1, then one can take s = 2, 

• if a = 2 and 2 < d < 3, then one can take s = 1, 

• if cr >3 then one can take s = a. 



As emphasized in [2], in some cases, the global existence of smooth solutions is already known 
for ([1]). For example, in the quintic case, u = 2, global existence is known for d < 3 (see [6] for 
the difficult critical case d = 3), so that only the bound (jlOp is interesting. Nevertheless, Theorem 
[T] may be also applied to cases where ([I]) is H^ super-critical {a > 3, d = 3 for example) and 
hence the fact that it is possible to construct a smooth solution on a time interval independent of 
e is already interesting. The main ingredient used in [2] is a subtle transformation of ([T]) into a 
perturbation of a quasilinear symmetric hyperbolic system with non smooth coefficients when a > 2. 



The first aim of this paper is to prove that the estimate (jlOp holds true for every s, every 
dimension d and every nonlinearity / which satisfies the following assumption: 



{A) 



/gC°°([0,+oo)), 



/(O) = 0, 



/'>0 on (0,+oo), 



3nGN*, /(")(0) /O 



Note that we allow /' to vanish at the origin. The assumption {A) takes into account in 
particular all the homogeneous polynomial nonlinearities f{p) = p" but also nonlinearities under 
the form f{p) = p"'^ + p""^ or j^ for example. Our result reads: 

Theorem 2 We assume (A), and consider an initial data ^ with Oq, (Pq in H°° such that, for 
some real-valued functions (c/3o,«o) £ H°^ , we have for every s, 



Cg - ao\ 



H" 



0{e) 



and 



V'o-V'o 



H^ 



0(e). 



Then, there exists T* > such that ([7]) with initial value (00,970) has a unique smooth maximal 
solution {a,ip) £ C{[0,T*[,H°° x H°°). Moreover, there exists T G (0,r*] such that for every 
e G (0, 1), the solution ^'^ to (HI)-© exists at least on [0,T] and satisfies for every s 



sup 

£6(0,1] 



^^ exp ( 03) 



< +00. 



L°°{[0,T],H^) 

More precisely, there exists ip^ = ip + Oh°°{£) such that, for every s, 



*'^exp ( c/5^) - a 



0(e) 



L°°{[0,T],H^) 



(11) 



Let us give a few comments on the statement of Theorem [2j 

At first, note that Theorem [2] contains a result of local existence of smooth solutions for ([9]) 
in the case of non necessarily homogeneous nonlinearities satisfying (A). Since (a, Vc/9) solves a 
compressible type Euler equation, the case of a homogeneous nonlinearity was studied in [15], and 
we thus give an extension of this result to smooth non-linearities satisfying assumption (A). A 
precise statement of our result with the required regularity of the initial data is given in Theorem 
H] below. The new difficulty when / is not homogeneous is that the nonlinear symmetrization does 
not seem to allow to transform the problem into a classical symmetric or symmetrizable hyperbolic 
system with smooth coefficients. 

The correction of order e that we have to add to the phase to get the estimate (llip is expected. 
Indeed, a perturbation of order e in the phase modifies the amplitude at the leading order. 

Our approach to prove Theorem [2] is completely different from the one of [2] and [9] . We do not 
work any more on the system ([6]) or any reformulation of ([T|) into a perturbation of a quasilinear 



symmetric hyperbolic system, but directly on the NLS equation ([T]). Basically, we first prove the 
linear stability for ([T]) in arbitrary Sobolev norms of highly oscillating solution of the form ae^^''^ 
and then use a fixed point argument to prove the nonlinear stability. The crucial estimate of linear 
stability of highly oscillating solution is given in Lemma [T] and Theorem [3l 

This actually allows to justify WKB expansions up to arbitrary orders (see Theorem [5]). Since 
we deal in this paper with sufficiently smooth and in particular bounded solutions, the assumption 
{A) can be replaced by a local version where we assume that /' > on (0, /3) with /5 independent 
of e if the initial datum verifies |aoP < /3- Indeed, since c^ takes it values in the (weak) hyperbolic 
region of the limit system ([7]), there still exists a local smooth solution of ([7|) defined on [0, T\ for 
some T > and the stability argument leading to Theorem [2] still holds. Consequently, our result 
can also be applied to nonlinearities like /(p) = p"^ — p^^ for every (T2 > o"! provided |aoP < /3 < 1. 
Note that when c72 is too large, the classical global existence result of weak solutions (see [8]) for 
([1]) is not valid and hence it does not seem possible to use the modulated energy method of [I] , [H] 
to investigate the semi-classical limit. 

Finally, the last advantage of our approach is that it can be easily generalized to the case of a 
domain with boundary and to non-zero condition at infinity. This will be the aim of the second 
part of the paper. We shall restrict ourself to a physical case, the Gross-Pitaevskii equation, i.e. 
f{p) = p — 1. The generalization to more general nonlinearities satisfying an assumption like {^A) 
is rather straightforward. This simplifying assumption is only made to avoid the multiplication of 
difficulties. Again to avoid too many technicalities, we restrict ourselves to the simplest domain 
= M^ = M'^-i X (0, +oo). For x G M'^, we shah use the notation x = (y, z), y E M'^-\ z > 0. We 
add to ([l]) the Neumann boundary condition 

a,^^(t,y,0)=0. (12) 

We also impose the following condition at infinity 

-it +i j, |x| ^ +00, (13) 

that we can write in hydrodynamical variables 

|^'^(t,a;)| ^1, u'^{t,x) ^ u°°, |x| ^ +00, 

where u°° is a constant vector. This condition appears naturally when we study a moving obstacle 
in the fluid. Indeed, if we start from ([1]) with the Neumann boundary condition on an obstacle 
moving at constant velocity and fluid at rest at inflnity, then we can use the Galilean invariance of 
([T|) to transform the problem into the study of ([1]) in a fixed domain but with the condition (J13p 
at infinity. 

This problem with such boundary conditions is physically meaningfull since it can be used to 
describe superfiuids past an obstacle (we refer to [16] for example). The semiclassical limit e tends 
to zero was already studied in [T3] by using the modulated energy method. The limit ([S]) was proven 
with (/9, u) the solution of the compressible Euler equation with boundary condition u ■ ti/qq = 0, n 
being the normal to the boundary. Note that the result of |14j is restricted to the two-dimensional 
case only in order to have a global solution in the energy space of ([T|) . By using more recent results 
on the Cauchy problem, [3], one can also get the result in the three-dimensional case at least when 
ii°° = 0. Our aim here is to give a more precise description of the convergence which takes into 



account boundary layers. More precisely, since the solution of the Euler system @ cannot match 
the Neumann boundary condition dza{t, y, 0) = 0, a boundary layer of weak amplitude £ and of 

size £ appears. They are formally described for example in p!6]. WKB expansions ^^ = a^e* ^ are 
thus to be seek under the form 

m m 

a^ = a^ + Y,e'{aHt,x)+A\t,y,^)), <^- = / + J]e'=(/(t, :c) + cl>'=(t,y, ^)) (14) 

fc=i fc=i 

where the profiles A^{t,y,Z), ^^{t,y, Z) are exponentially decreasing in the Z variable and are 
chosen such that 

d,a''{t,y,0) + dzA^+\t,y,0) = 0, d,v''{t,y,0) + dz^''+Ht,y,0) = 

so that the approximate WKB expansion xj/"^"^-" = a^ exp (-93^) matches the Neumann boundary 
condition ()12p . Our result (Theorem[6]) is that under suitable assumptions on the initial conditions, 
we have the nonlinear stability of WKB expansions: in particular we have the existence of a smooth 
solution for ([1]), (J12p . (J13p on a time interval independent of e and the estimate 

ii^'e-^^-a^||^,,^<e. (15) 

Note that it is necessary to incorporate the boundary layer £A^ in order to get (jlSp since its gra- 
dient has amplitude one in L°°. The case of Dirichlet boundary condition which is also physically 
meaningfull, we again refer to [TB], seems more complicated to handle as often in boundary layer 
theory in fluid mechanics since the boundary layers involved have amplitude one. This is left for 
future work. 

The paper is organized as follows. In section [21 we prove the linear stability in H^ of an 
approximate WKB solution of ([T]) under the form a^ exp (i— ) in the case Q, = M . This is the 
crucial part towards the proof of Theorem [2j Next in section [3l we give the construction of a 
WKB expansion up to arbitrary order and give the proof of the local existence of smooth solution 
for the compressible Euler equation with a pressure law satisfying (A) . In section HI we give the 
justification of WKB expansions at every order and recover Theorem [2] as a particular case. This 
part uses in a classical way the linear stability result and a fixed point argument. Finally, in section 
[5l we study the problem in the half-space with Neumann boundary condition. 

2 Linear Stability 

In this section, we consider a smooth WKB approximate solution ^" = a^ exp (i^) of ([T]) such 



that 



NLSi^"") = R^ exp [i^) , (16) 



where 

NLS{^) = i£dt^ + — Al- - ^fil^l"^). 



Moreover, we also set 

R^^dtip' + l\V^'\' + f{\a'\'), (17) 

Ra^dta' + V^' ■Va' + -a'Aip', (18) 

so that 

R^ = -d'R^ + ieRa- 

Looking for an exact solution of ([1]) under the form 

If ^ = ^« + u; e'^ = (a^ + w)e'^, 
we find that w solves the nonlinear Schrodinger equation 

/ I \ e^ 

i£{dtw + u^ -Vw + -wV ■ u''] +—/\w-2{w,a!')f'{\a''\^)a!' = R^w - R"" + Q{w), (19) 

where we have set 

and the nonlinear term Q{w) is defined by 

Q(^) ^ (a- + ^) (/(|a- + ^|2) _ /(|„.|2)^ _ 2(^, a^)/'(|a^|2)a^ (20) 

Of course, i?^ will be very small and Rip (and Ra) are to be thought small (at least 0(e)) for 
applications to nonlinear stability results. Nevertheless, in this section the exact form of these 
terms is not important. The way to construct an accurate WKB solution will be explained in the 
next section. 

Remark 1 If we work with a non-linearity / such that /(|ylp) = 0, we can impose a non-zero 
condition at infinity such as oq G ^ + H°° and Vipo G U°° + H°° . Since we can still look for the 
perturbation w in H^, this does not affect the proofs. 

Since we expect the correction term w to be small, we shall only consider in this section the 

linearized equation 

dw 

where the linear operator C^ is defined as 



ie'^+C''w = RpW + F'', x G M*^, (21) 



2 

C^{w) = —/\w + ieu^-Vw + -wV-u^-2f{\a^\'^){w,a^)a^ 

and (•, •) stands for the real scalar product in C ~ M^. In this section, F^ is considered as a given 
source term. Of course, for the proof of Theorem [21 we shall apply the result of this section to 

F' = -R' + Q{w). (22) 



We notice that the first and last terms in C^w are formally self-adjoint, and that the quadratic 
form (in H^) associated to the operator 



IS 



/ 



f2 






d_ 
di 



It is then natural to consider the (squared) norm / (^'w,S^{w)^ as a good energy for the linearized 
equation (j2ip . Consequently, we introduce the weighted norm 

N'{w) = - [ e2|Vu;|2 + 4/'(|a"|2)(u;,a^)2 + Ke2|^|2 

for every K > (K will be chosen sufficiently large only in the next subsection). 

Our first result of this section is a linear stability result in the energy norm N^{w). 

Lemma 1 Assume that u'^ : [0, T] x M'^ ^ W^ and a^ : [0, T]xW^ ^ C are smooth and such that 

Let w G C^([0,T], i/2) he a solution of (1211) . Then, there exists Cm depending only on d, f and M 
such that for every e £ (0, 1], the solution of ()2ip satisfies the energy estimate 

N'{w{t)) < CmU + ^||i?a(t)|Loo + ^\\RAt)\\w^,o. + ^||^^(OlLoo V"(«^(i)) (23) 
+ ll^'W|li2- / -f'{\a'\'){w,an{a',iF^)+ [ {eAw,iFn. 

Note that it is very easy to get from (j23p and the Gronwall inequality a classical estimate of 
linear stability. Indeed, assuming that Ra = C'loo([o,t],l°°)(^) and Rip = 0^oo([o,T],iyi'°°)(£^) (which 
is true if {a^,ip^) come from the WKB method), we infer from a crude estimate for the two last 
terms in ^ that for < t < T, 

J^N^Hi)) < CN%wit)) + ^\\F%t)\\l„ 
which gives for < t <T 

N%w{t)) < e^*(iV^(^(0)) + 1|* \\F%r)\\%, dr), 

which is a more classical result of linear stability in the energy norm N'^{w) since the amplification 
rate C is independent of e. Nevertheless, to get H'^ estimates and the best nonlinear results as 
possible, it is important to have the special structure of the two last terms in (j23p . 

Modulated linearized functionals like N^ were also used in asymptotic problems in fluid me- 
chanics, see |10j for example. 



2.1 Proof of Lemma [T] 

The norms L°°, W^'°°, L^ ... always stand for the norms in the x variable. 
At first, since S^ is self adjoint, we have 

^f {S'w,w)=[ 2{S'w,dtw)+2dt[f'{\a'\^)]{w,a'f + Af'{\a'\^){w,a'){w,dta'). (24) 

Next, we use ([2T]l to express dtw as 



e ^ 2 ^ e e 



to get 



2 / {S^w,dtw) =2 / (^—^w-2f{\a'\'^){w,a^)a',u'-Vw + -wV-u^ + -R^w + -F'^). (25) 



We shall now estimate the various terms in the right-hand side of (j25p . Integrating by parts, we 
get 






^^" ^ 'eHVwf.^+e'^Wwf 



< e\\VRip\\Loo\\w\\i2\\Vw\\]^2 

< -\\R^\\wi,ooN%w). 

Note that we have used that R^p is real-valued and thus that 

{Vw,iR^\/w) =0 
for the first equality. We also easily obtain by integration by parts that 

< CmN'{w). 

In the proof, Cm is a harmless number which changes from line to line and which depends only on 
M. In particular, it is independent of e. Moreover, we can also write for k = 1,- ■ ■ , d, 

f {dlkW,u'-Vw) = - f n^.V^^- f{dkw,dku'-^w) 

jR'i jR'i ^ J 

and hence, we immediately infer 

/ (e^Aw, u^ ■ Vw] < CmN^{w). 



Furthermore, from the inequaUty 2ab < a? + b'^^ there holds 



1 



< ^\\RM- f f'{\a'\'){w,an'+e^\w\^ 

< ^WRM-'N'H- (26) 



Consequently, we can replace ([25]) in (p^ and use the above estimates to get 

^ [ {S'w,w) = [ Af'{\a'\'^){w,a')({w,dta')-{u'-Vw + -wV-u^,a^)) (27) 
+2/ dt[f'{\a'\')]{w,a'f + Eu 
where -Ei satisfies the estimate 



El < Cm{1 + -\\R4^i^o. + ^\\R4lo.]N%w) (28) 



e "" ' e^ 



u' 



-- / /'(|a1')(«i,a^)(a^in+ / {eAw,iF'). 
To estimate the first integral in the right hand side of (f27l) . we use the equation (fTSJ) to get 
4 / f'{\a'\^){w,a')({w,dta')-{u'-Vw + lwV-u',a')) 

jR'i ^ 2 / 

= 4 / f'{\a'\^){w,a')({w,Ra) - u' ■ V{w,a') - {w,a')V ■ u') 

= 4/ f{\a'\')iw,a'){w,Ra)-2 f f'{\a^\'')u^.V{{w,a'f)-^f n\a'\''){w,a'fV 

jRd JM.'l JR'i 

= 4/ f{\a'\'){w,a'){w,Ra) + 2 [ {w,a^fu'-V[f{\a'\')]-2[ /' (|a^|2)(^,a-)2 V • n^ 

JRd JR'* JR<^ 

To get the last line, we have integrated by parts the second integral. Note that the last term is 

bounded by CmN^{w), and, as for (|26p that the first integral is bounded by \\Ra\\L°°N'^{'w)- 

Consequently, we can replace the above identity in ([271) to get 

^ / {S'w,w)= [ 2iw,a'f(dt + u'-v)f'{\a'\^) + Ei+E2=:I + Ei + E2, (29) 

at Jud jRd \ / 

where £'2 is such that 

E2<CM(l + ^\\Ra\\Lo?jN%w). (30) 

To estimate /, we use again the equation (fT8|) which gives 

(a* + u' • v)/'(|a^|2) = 2r(|a^|2)(a^ dta' + u' • Va') = 2f"{\a'\^)(^Ra -^a'V-u',c 

10 



and hence we find 

I<C f \a'nf"{\a'\')\{w,a'f + 4 I \a'\\f"{\a'\')\{w,an^\Ra\. 



To conclude, we shall use the assumption (A). By defining n G N* the first integer such that 
y(")(0) ^ 0, we see from Taylor expansion that 



f'ip) = P^'\{P) 
for some smooth positive function q on [0, +00). In particular, since g > 0, we have 



(31) 



Pf (p) 1 , Qip) ^ roo(\n , ^\ 

f [P) Q[P) 



which implies 
This yields 



\pf"ip)\<CMnp) 



for 0< p<M. 



(32) 



\a'\^\f"{\a'\^)\{w,a'f<CM / {w,an^f'{\a'\^)<CMN^w), 
where, again. Cm depends only on M. In a similar way, we also obtain 

{w,an^\a'\ \f"i\a'\')\ \Ra\ < ||^a|L^ / H ■ \{w,a')\ ■ \a'\^ \f"{\a'\'' 

Cm 



< 



\Ra\ 



[e\w 



I) 



>,a-)/fG^ 



< ^l^-IL^^^M 



Consequently, we have proven that 



I<CM{l + -\\RaL^)N'{w) 



(33) 

To get the result of Lemma [TJ it remains to perform the L^ estimate. Taking the L^ scalar product 
of ([^T|) with iw and using that 



[w,u 



we get 



d /e^ .. ..0 



Vw + -wV -u^) = -V ■ {Iwfu"), 



e{F^,iw) + 2e [ f'{\a'\'^){w,a^){a^,iw). 



Note that we have once again used that R^p is real-valued and hence that {R^w, iw) = 0. The first 
integral is clearly bounded by N^{w) + |-^^|| r2 whereas for the second one, we have 

-2ef'{\a'^\^){w,a'){a',iw)<CM [ (n\a^\^){w,a^f + e^\w\^) <CMN'{w{t)). 

As a consequence, we get 

d /e^ 



:^(Tll^lliO^^^^^'H + ll^1li- 



dt\2 
Finally, we can collect ([28]), ([29]), ([30]), ([MD and dHl) to get ([23]). This completes the proof. D 



(34) 
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2.2 Higher order estimates 

Since our final aim is to prove Theorem [2] by a fixed point argument, we also need to have 
H^ estimates for s sufficiently large for the solution of the linear equation ()2ip . This is the aim 
of the following. Note that the term — 2(i(;, a'')/'(|a^p)a^ in (J19p can be seen as a singular term 
with variable coefficients. Consequently, a crude way to get H'^ estimates is to apply el"l5" to 
the equation, the weight el°l being used to compensate the singular commutator when we take the 
derivative of (J19p . and then to apply Lemma[T]to the resulting equation. Nevertheless, it is possible 
to avoid the loss of e'°' with more work by using more clever higher order modulated functionals. 
If s € N, s > 2, we define the following weighted norm, where a G N*^ are multi-indices 



Nl{w) = Y^ N^{d''w)+K\\R(iw\]j.-2 (35) 

\a\<s-l 



\a\<s-l 

In this section, we shall use that 



a^ = a^ + ea' 



with a^ real- valued and 



Note that this allows to write 



sup ||o''||l°=([o,t].vk='°°) < C. 

eG(0,l] 



and hence by choosing K sufficiently large {K > C) we get the lower bound 

Nl{w)>\ Yl ^'i9»+ E / n\a'\'){a'f\Re d-w\' dx (36) 



a|<s— 1 |q|<s— 1 

Note that we also have the equivalence of norms: 



wY^, < -Nliw), Nliw) < Ci\a'\w^-i,^)\\w\\'^s + ||Re w\\j,..,. (37) 



e 
The main result of this section is : 

Theorem 3 Let < T < oo, s £ W , f satisfying (A) and w G C'^{[0,T],H^) a solution of dH 
with u" : [0, T] X M'^ ^ R'^ and a^ : [0, T]xW^ ^C such that 

M = sup I ||U ||^oo([o,T],VF=+i'°°(R'')) ~'~ in llL°=([0,T],VF='°°''rod~i^ ) ^ +00. 



Assume finally that, for some a^ G L°°([0, T], 1^*'°°(M'^)) real-valued, a^ writes 

a^ = a° + 01^^,00(6) (38) 

uniformly on [0, T] . Then, there exists C , depending only on d, f and M , such that 
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Remark 2 In view of (j38p . a^ is real up to 0(e), hence, in the integral in the right-hand side of 
()23|) . the real and imaginary parts of F^ do not play the same role. This explains that the estimate 
is better for Re-F^ than for ImF^. As a matter of fact, for s = 1, Theorem [3] follows immediately 
from Lemma [Hand (jSSh . 

2.3 Proof of Theorem H 

We estimate separately the two terms in Ng{w), when s > 2 (otherwise, the result follows from 
Lemma [Has we have seen). Let us set 

Y,{w) = ||Re w\\'jjs-2- 

Note that we have 

S(u;) < Nliw). (39) 

In the proof, C is a constant depending only on d, f and M. 
We shall first prove that 



|S(«;) < C(l + }^\\R^\\^^,,^^)n!{w) + C\\F^\\l..^, + f 



-E{w) < C(l + -\\R^\\^^^,^\nI{w) + C||F^||^,_, + -||lmF^||^,_,. (40) 



For a G N , we have 

IF 7 7 

dt{d"w) + u' ■ V(9"u;) = -A{d'^w) - -d'^F' - -d''{R^w) (41) 

^ Co 

--a"(/'(|a"|2)(a^u;)a^)- [a",n" • V]u; - ^ 5°('u;V • n"). 
Next, by taking the real part of (j4ip . we get 

at(a"Reti;) +u^-V(a"Retf) = - [9", u" • VJRe tu - ^ ^"(Re u; V • n^) + 7^^ 
where 

7^^ = Re(-A(a"w;) - -^"F^ - -^"(iJ^u') - -a"(/'(|a^|2)(a^u;)a=)) . (42) 



By using (|3S|) . we have 

Im 9^a^ = 0(e), V7, I7I < \a\ 

and 

|(5'3a^a^^t;)| < C|Re d^'w\+e\d^w\ (43) 

for every /?, 7. Consequently, we immediately obtain for every a, \a\ < s — 2, 



11 _ 11 / 11 11 -tXjin TA/s — 1,00 II II II -p.. II II II \ II T 7— ?p II 

f/v L2 < C 5hf^Lf/s H ^- — - — ■ — hf^ f7s-2 + Re w\us-2 + ^w? 775-2 H — Im F \tj. 



£"■ J £ 



3-2 



< c(l + "-^^"y-^- )Ar|(^)l + l||ini F^ll^.-.. 
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Consequently, the standard L^ energy estimate for (|i2]) gives 



Note that we have used that 



^||Re 5-||i. < C(l + ""^"^;-"°° )iV|H + ^||Im F^||2^._.. (44) 



/ (u' ■ V(a"Re w),d'^Re v?j = -]- /(V • n^)|a"Re w\ 
Consequently, (00]) is proven. 



The next step is to estimate N^{d°'w) for |q| < s — 1. By applying d'^ to (I2ip . we get 



ze 



+ £' {d'^w) =R^d''w + F', (45) 



where 
with 



dt 

C° = 2 9°(/'(|a"|2)a^(Tz;,a")) -2/'(|a"|2)(9°^z;,o")a^ 
P" = -ie[a",u^-V]u;--[a", V-u%. 

To estimate N^{d°'w), we shall use Lemma [TJ Towards this, we need to estimate the commutators 
in the right hand side of (I45p . For |a| < s — 1, the following estimates hold for C° and D": 

\\[d'',RM\m ^ C\Mw^^4M\l^ < ^\\RAw.,ooN!iw), (46) 

P%r <Ce''\\w\\l,<CNliw), (47) 

||(i/'(|a^|2)la^P")||^, <Ce2Ar|(«;), (48) 

\\C%, <CN!{w), (49) 

||(m^C")||^, <Ce'N!{w). (50) 

The estimates ()46p and (j47p follow easily from (j37p . For (j48p . we note that 
-(m",P") = -(a^[a",u"•V]u;)-^(a^[9°,V•n^]u;) 

since u"" is real. Next, we can use ()38p and (j43p again. In particular, in the above expansion, the 
terms (a^, d'^w) are bounded in L^ by S(u') + e^||t(;||^3_2 and thus by N^{w). Similarly, the terms 
(^a'^,Vd'^w) are bounded in L^ by A^f(w) if I7I < s — 3. Consequently, we get 

||(^/'(|a^|2)la^P°)||^, <C( J^ / /'(|a^|2)(9V a^)' + iV|(t«)) < CiV|(^) 
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which yields (j48p . Next, we turn to C". The Leibnitz formula gives 

C"= Yl *d^[f'{\a'f)]{d"w,d'^a')d''a', (51) 



a < a, 

a + f3 + \ + fi = a 



where * is a real coefficient depending only on d, (5, A and /i. Since \d\ < \a\ — 1 < s — 2, we can 
use again ([38]) through (|i3|) to get that 



||C"||i2 < C[i:{w)+e'\\w\\ijs'j < CNliw). 

Since (ia'^,5^a^) = 0{e) thanks to (j38p . we also get (j50p . For the H^ norm, the same argument 
yields 

\\C''fHr<c(^{w)+e'^\\wfHs+ V [ \d^[f'{\a'\'^)]{d^w,d'^a')d^'a' 

^ II-, J^"^ 

7| = s — 1, 

l/3 + A + ^l = 1 

To estimate the last sum, we first consider the terms with /3 = 0. They are always bounded by 

C f [f'i\a^\') + \a^'\n\a^')\]{d'^,af 
with I7I = s — 1 and hence, thanks to ([32|) . they are bounded by 

C [ fi\a^\'){d'^w,a'f 

jRd 

and hence by Ng{w). Next, we consider the terms with |/?| = 1. Since A = /i = 0, we have to 
estimate terms like 

Ti= [ f{\a^\'){d^w,df'af\a^\'. 



By using again (|38|) and (jl3|) . we get 

Ti<C f /'(|a"|2)|a°p|Re a^w;^ + Ce2||y;||2^^_^ 

and hence, by using (j36p . we finally obain 

Ti < CiV|H. 

Consequently, (|49p is proven. This ends the estimates of the commutators. 
We are now able to establish: 

j^N'{d''w) <C {1 + l||i?^||^^_,_^ + i II i?, II ^^)iV|H 
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Indeed, from Lemma [H we deduce 



dt 



+ II^1IL + - / f'{\a'\^){d''w,a%ia',F')- f {zeAd^w^F^, (53) 



To estimate the right-hand side of ()53p . we first estimate ||F^|L2- Combining ()46p and ()47p with 
^9l). we infer 

II^IIl^ < ll^1li.-i +^(l + ^ll^^ll'H/.-i,o.)iV.^H. (54) 

Next, we turn to the term 



which sphts as four integrals. For the first one, by (I50p and Cauchy-Schwarz: 

/'(|a^|2)(9"u;,a^)(ia^C") <c( f /'(|a^|2)(5"^,a^)2) ^ iVf (i.)l < CNHw). 
For the second one, we use (|48p and Cauchy-Schwarz, which gives 

/'(|a^|2)5(9°u;,a^)(i/'(|a^|2)5a^P") < CNI{w). 



4 

£ 

For the third integral, we simply write, using once again (j38p 

-\\{za',d-F')\\^,<C\\F'\\^^^,+-\\lmF'\\^,_„ 
which yields by Cauchy-Schwarz 



/ f'{\a'\^){d''w,a'){ia\d''F')<CN!{w) + C\\F'\\l^_, + ^\\lmF'\\l,_,. 



Finally, for the fourth integral, we have by ([36 



- / f'{\a'\^){d''w,a%ia\[d'',RM<-\Mw^^^-^^"H- 
By summing these estimates, we find 
- / /'(|a^P)(5"u;,a^)(ia^F0 <c(l + -||ii^||^._,,^)iV|(u;)+C||F^||^._, + -^||lmF^||^,_,. 

£ J]^d \ £ / £ 

(55) 
Finally, we handle the term 

- / {i£Ad''w, F^) = - I (ieA9"u;, C" + P" + ^"F" + [9", R^]w) . 
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By using an integration by parts, we have 

- / {ieAd''w,F') < ||C"||' + ||P"||' + \\[d'',R^]w\\l, + \\F'\\%^ + CNl{w) 

thanks to (gT]), (06]) and (09]). 

Consequently, we can collect the last estimate and ([M|) . ([55]) . ([53]) to get ([52|) . This ends the 
proof of Theorem [S] 

3 Construction of WKB expansions 

In this section, we construct an approximate solution of ([T]) using a WKB expansion. The first 
step is to prove the local existence of smooth solutions of the limit hydrodynamical system. 

3.1 Well-posedness of the limit system 

We consider the system 

dta + u ■ Va H — aV ■ u = 

2 (56) 

, dtu + u-Vu + V{f{a^)) =0, 

which is only weakly hyperbolic, with the pressure law / satisfying assumption (A), and the initial 
condition (a,u)|t=o = {a^,u^)- 

Theorem 4 Asswne that f satisfies assumption (A), and denote 6y n € N* the first integer such 
that /(")(0) / 0- Then, for s > 2 + d/2, for every (oq, uo) e ^^"^ x H" with ao G M and a[J E H" , 
there exists T > and a unique solution {a,u) of (I56p such that {a,a'^,u) G C{[0,T], H^^^ x H^ x 
H')nC^{[0,T],H'~^ X H'-^ X H'-'^). 

Let us remark that if n = 1, then /'(O) > and thus /' > in [0, +oo) (by {A)). In this case, 
([56l) is symmetrizable (using the symmetrizer S = diag(l, ^/'(a^), ..., jf'{a^)) used in [9j) and the 
existence and uniqueness for ([56]) follows easily. 

Proof of Theorem [4] 

The first step is to rewrite the system by using more convenient unknowns. At first, we notice that 
thanks to {A), we can write / under the form 

/(p) = P"/(P), 

with / smooth on [0, +cxd) and such that /(O) 7^ 0. Next, since we have by assumption /(O) = 
and f'{p) > for p 7^ 0, we also have that f{p) > for p > 0. This implies that f{p) > for p > 0. 
This allows to define a smooth function /i on M by 



1 



h(a)^a[/V)]^"- (57) 
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Note that h(a) 7^ for a 7^ 0. It is usefull to notice that we can also write h under the form 

h(a) =sgn(a)/(a^)2^ 

and hence that we have 

h(a)2- = /(a2). 

Furthermore, since /' > and /(O) > 0, we deduce that h'(a) > for a 7^ and that h'(0) = 

1 
[/(O)] ^" > 0, so that h' > on M. Thus h is a smooth diffeomorphism from M to h(]R). In 

particular, this allows to define a smooth positive function c on h(]R) such that 

-ah'(a) =h(a)c(h(a)), \faeM. 

With this definition, {h,u), with h = h(a), solves the system 

dth + u-Vh + hc{h) V • u = 

(58) 
^ atn + u- Vu + V(/i2") =0. 

Since a is in W if and only if h is in H^ ^ we shall prove local existence of smooth solution for the 
weakly hyperbolic system (j58p . As we shall see below, the nonlinear symmetrization method of 
|15| does not allow to reduce ()58p to a symmetric or symmetrizable system with smooth coefficients 
except in the case that c{h) = c{h"') with c smooth . Nevertheless, it will be still possible to use the 
same idea to prove the existence of an energy estimate with loss for the system ()58p . When we are in 
such a situation, the simplest way to construct a solution is to use the vanishing viscosity method. 
Indeed, this approximation method allows to preserve the nonlinear energy estimate verified by 
(f58l) . We thus consider for e > the system 

dthe + u^ ■ V/ie + hf:c{he)V ■ u^ = e Ahf: 

(59) 
dtU^ + U,- VUe + V (/i^") = e AUe . 

The local existence of smooth solutions for this parabolic system is very easy to obtain. Moreover, 
we note that h^ remains nonnegative if the initial datum /i|f=o is nonnegative. In the following, we 
shall only prove an H^ energy estimate independent of e for this system which ensures that the 
solution remains smooth on an interval of time independent of e. The final step which consists in 
using the uniform bounds to pass to the limit when e goes to zero to get a solution of ()58p is very 
classical and hence will not be detailled. In the proof of the energy estimates, we shall omit the 
subscript e for notational convenience. 

As in the work of [15], we introduce the unknown H = K" = a^f^a^)^. Note that by definition 
of /i, H is in H" as soon as a" is in H'^ . We get for [H, u) the system 

dtH + U-VH + nHc{h) V ■u = en/i"-^A/i = e (aH - n{n - l)h'^-'^\Vh\'^) 

^ ' (60) 

dtu +u-Vu + 2HVH = e An. 

Note that it does not seem possible to get a classical hyperbolic symmetric system (in the case 
e = 0) involving only H and u as in the case of homogeneous pressure laws considered in |15] . 
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Indeed, the coefficient c(/i) = c{H~) is not (in general) a smooth function of H. Nevertheless, it 
will be possible to prove that the system with unknowns (/i, H, u) though only weakly hyperbolic 
(when e = 0) satisfies an energy estimate. We notice that the symmetrizer 

5^diag(l,^c(/i)/rf), 

which is positive since c{h) is positive, symmetrizes the first order part of (j60p . We shall first 
perform an H^ energy estimate (s > 2 + d/2) on (I60p but we have to track carefully the dependence 
on h in the energy estimates. 

To prove our H^ energy estimate, we shall make an extensive use of the following classical (see 
|18j for example) tame estimates 



Ifalijk < Ck[\\f\\Loo\\g\\fjk + ||/||j:^fc||fi'||iooj, (61) 

||9"(/5) - /5"5|L2 < Ck[\\f\\j,,\\g\\^^ + ||V/||^^||5||^._i), |a| < k, (62) 

\\F{u)\\j^, < C{\\u\\^^){l + \\u\\^,) (63) 

if F is smooth and such that F{0) = 0. 

At first, we notice that {d'^H,d"u) for |a| < s solves the system 

' dtd'^H + u-Vd''H + nc{h){V ■u)d''H = e (a^"/? - n(n - l)a"(/i"~2|V/i|2)) 

-[d",u\ • VH - n[a", Hc(h)]V • u 

dtd^'u +u-va"u + 2HVd°'H = eAd'^u - [d'',u] -Vu- [d'^,2H]VH. 

By using (j62l) to estimate in L^ the commutators in the right hand-side, we get in a classical way 
by integration by parts 

^\l f |a"i7|2 + ^c(/i)|5"u|2j +e / |Va"//|2 + ^c(/i)|Va"n|2 (64) 

dt 12 J^d 2 J J^d 2 



<Co(||(/l,n)||^,,^)||y|||,.+C" + eP" + 7^^ 
where V = {H,u), Cq is a non-decreasing function depending only on /, s and d, and 
C" = -n / (a"F) [d'',Hc{h)]{V ■ u), 

P" = -j/ c'(/i)((V/i-V)9"n)-a°n-n(n-l) / a"(/i"-2|V/ip) 5"//, 

7^" = ^ /" c'(/i)5i/i|a"up. 

We have singled out the three terms above since they are the ones involving h which must be 
estimated with care. Note that the estimate of C" will be crucial since this term involves high order 
derivatives of h. Next, we can integrate ()64p in time, sum the estimates for |a| < s and use that 
c{h) > 0, hence nc{h)/2 > ^ . ii A \ to obtain 

\\V{t)\ls+ef\VV{T)\lsdT (65) 

JO 

< C'i(||/i||L-)(||y(0)|||,. + jJ|*C'o(||(/i,n)(T)||t^i,^)||y(r)||f,. +C(t) + eP(T) +7^(T)dr) 
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with 

q|<s IckI^'^ |ci|<'5 

Estimate for C. We claim that 



C<Co(||(/i,u)||,^i,. 



H" + ll"-llH-s-i 



II?.. 



(66) 



The crucial point is that this estimate only involves the H^ ^ norm of h. This will allow to conclude 
by using that for the first equation in (j60p . the H'^~^ norm of h is controlled by the H'^ norm of u. 



By using the commutator estimate (I62p . we have 

C < C\\H\\Hs(\\Hc{h)\\H^\\V-u\\Lo^ + \\V{Hc{h))\\L--\\V-u\\Hs-i) 
< Co{\\{h,u)\\wi,oo){\\VfHs + \\H\\H^\\Hc{h)\\Hs). 

To estimate the last term, we use that H = h^, which yields hdiH = nHdih, thus 

1 



di(Hc{h)j = c{h)d^H + c{h)Hdih = c{h)diH + -c'{h)hdiH. 
Consequently, by ([6T]) . ([63|) . we get 

\\Hc{h)\\H^ <C\\c{h)VH\\Hs-i+C\\c'{h)hVH\\Hs-i <Co{\\ih,u)\\wi.oo)(lH\\H^ + \\h\\Hs^^y 
and (1661) follows. 



Estimate for T>. The term D involves derivatives of u of order < s + 1, and we shall use the 
energy dissipation in ([M|). We prove that 



Ci 



L-)eP< -ejVVfHs +e 



lyi 



h + Wh\\l. 



(67) 



We have, on the one hand. 



[ c'{h)Vh-Vd'^u-d'^u 

JRd 



<Cr 



0{\\n\mn,ac^ 



Vli||f/s WuWh" < Cr 



oUKi|H4/i,oc, 



WvvWhsWv 



H" \\y WH"- 



On the other hand, for the second term (which vanishes if n = 1), after one integration by parts 
when |a| > 0, we get 

n(n-l) [ a"(/i"-2|V/i|2)a"// < C\\VH\\Hs\\h''-^\Vh\'^\\Hs-i 

< Coi\\h\\wi.o.)\\VH\\H4^h\\Hs-i, 



and if a = 0, since H = h^ and s > 1, 

n(n-l) / h^'^\Vh\'^H 



n — 1 

n jRd 



|Vif |2 < CWHWJjs 
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Consequently, 

eV < eCo(||/i||^i,oo) \\VV\\hs {\\V\\h^ + ||V/i||h-i) + e C\\V\]js, 
and (j67p follows from the standard inequality, for a, b, 6 > 0, ab < 9a^ + ^. 
Estimate for TZ. We prove that 

Ci{\h\\L^)n<\e\\VVfHs+CMh,u)\\w^,o.)\\V\ls. 



By using the first equation in (j59p for /i and an integration by parts, we find, as for the first 
term in 2?, 



9 n 



ia„.|2 



c'(/i)A/i|a°n 



4 ./rod 4 



2 ina„,|2 



c"(5)|v/iha"u 



7^° < Co(||(/l,u)||iyi,oo: 

< C7o(||(/i,u)||h^i,oo: 

< C7o(||(/i,u)||h^i,oo) (^||y||^. +e||vyb^ ||v 11^. 
Then, ([67|) follows as above from the inequality ah < 9a? + |g. 

Summing ([66]) . ([67]) and ([68]) . inserting this into ([65|) and simplifying by e ||Vy|||^s, we infer 



\v{t)\\s<c,{\h{t)u^)(\\vm 



Ms 



+ / C7o(||(/i,n)(r)||^i,oc) ||F(r)||l,. + ||Mr)||^.., + efV/iO 



II ff- 



dr . 



(69) 



To close the estimate, it remains to estimate ||/i||^s_i and e Jq || V/i||^s_i . We use the standard H^ ^ 
estimate for the convection diffusion equation (j59p which yields, as for (j64p . for |a| < s — 1, 



d 
di 



,aL|2 



\d'^h\ 



+ e 



\d"hf<Co{\\{h,u)\\^,^^){l 



\\H= 



+ \\h\\Hs-i\\u\\H'' • 



Summing for |q;| < s — 1 and integrating in time, this yields 



(70) 
Finally, we can combine ([69]) and ([70]) . to get 

\\v{t)\\h+\\Ht)\\%-r 



<Co(||(/i,n)||^^([o^^j^^,,^))(||F(0)||f,. +||M0)||2,._, + / ||F( 



(r)||l,. + ||Mr^"^ 



Hs-i dTJ. 

(71) 



Since i?^ ^ is embedded in T^^'°° for s > 2 + (i/2, we easily get by classical continuation arguments 
and the Gronwall lemma that the solution of (j59p is defined on an interval of time [0, T) independent 
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of e. Finally, (j7ip provides a uniform bound for {h, H, u) in H^^^ x i?* x /:f *, which allows to prove 
in a classical way that (h^,Ue) converges towards a solution of ()58p . This ends the proof of the 
existence of solution. 

To prove the uniqueness, it suffices to use the same method as above and perform an L^ energy 
estimate on the system satisfied by hi — /i2, ui — U2, Hi — H2. This is left to the reader. 

3.2 WKB expansions 

We now turn to the construction of WKB expansions up to arbitrary order. Let us first notice 
that in Theorem IU if the initial datum (ao,Mo) is in H°° x i^°°, then the solution {a^u) is in 
C°([0,r],iJ'*-i X H") for every s > 2 + d/2, with T independent of s > 2 + d/2. In other words, 
the existence time of the maximal solution in H°° x H°° is positive. This fact follows easily from 
([71]) and the Gronwall inequality (since H'^^^ C W^''^). 

Lemma 2 Consider ^q = aQC^^ol^ with Og G H'^ , tp^ G H'^ and that for some m G N, there exists 
an expansion 

m m 

a^ = ^e%^ + e™+iag, ^g = ^ e Vo' + e'"^ V^ (72) 

fc=0 fc=0 

with Qq G M, satisfying, for every s, 



£6(0,1) ^ ^ 

Let us denote T* is the maximal time of existence of a smooth (i.e. H°^ x H'^ ) solution {a^,ip^) 
for (j56p with the initial condition (aQjC^p). Then, there exists an approximate smooth solution of 
([1]) on [0,T*) under the form ^'' = a^e*'^"''^, with a^,(p^ G H^ and a^ complex-valued, solving 

(74) 
^ — + (Vv.^) • Va^ + -Av9^ - JeAa^ = ii^, 

wii/i i/ie initial condition (a^, (p^) , = (ag, (/jg) , and where, for every s and < T < T* , 

sup (\\K\\hs + II^^IIhO ^ C',,Te™+^ (75) 

[o,r] ^ ^ 

Fina%, forO<T< T* , a' verifies §Si): a^ - a^ = 0{e) in L°°{[0,T],W'''^). 
Note that ^" is indeed an approximate solution of ([1]) since 

ie^ + — A^'^ - *V(I^1^) = ( - ^e«^ + a'R"^) exp (i-^). 



By using the notation of section [2l we have R^ = —ieK^ + a'^R"}, hence 



sup||i?^||^, <C,e"+2, (76) 

[0,T] 
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Proof. 

As in [9j, we look for expansions 



Y^e'^a^ + e^'+^a 



^= = ^e*^/ + e"^+V 



m+l 



fc=0 



fc=0 



This yields that (a , i^ ) solves the nonlinear system 

which is just ([9]), and that for 1 < k < m, (a , (/? ) solves the linear system 



dt 



da 



+ 2fi\a^\'){a", a'^) + V^" ■ V^^ = S' 



^ + (V/) • Va'^ + Va° • V/ + ^ A/ + ^ A/ = 5^ 



9t 



(77) 



(78) 



where the source terms (5^,5^) depend only on (a-', v3-')o<j<fc-i, and S^ is complex-valued. 

We first solve ([77|) (that is ([7])) with the initial condition v?/j^q = v?o, O/j^g ~ ^o- -^y introducing 
ti'^ = Vt/j*^ and by taking the gradient of the first equation of ([771) . we find 



9taO + txO-VaO + -V-nO = 



,0^2 



(79) 



^ atn° + u°-Vn° + V(/(( 



0, 



which is the compressible Euler type equation considered in the previous section. By using Theorem 
m we get the existence of a smooth solution (a^,u^) £ H^^^ x H^ for every s on [0,T*) (with T* 
independent of s), with a^ real- valued. Finally, to get ip^, we can use the same argument as in [^. 
This yields 

^\t,x) = v.°(x) - I* (/((a0)2) + 1 |t.0|2)(r,x)dr. 

We now turn to the resolution of (j78p . We solve it with the initial condition [^ip^,a^^,„ = 
(^Jq, Oq). By introducing again u'^ = V(p^, we can take the gradient in the first line of ()78p to get 



_ ^iu'^ + n° • Vu'' + V(a°, /'((a°)2)a^) + n'= • Vu° = VS^. 



(80) 
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Again, since f'{{aP)'^^ can vanish, the symmetrization of this hnear hyperbohc system requires 
some care. We thus set 

j,k(, .j V2{n{a^f)f^a^ if n is odd 

^'"^^"1 V2a\f^Y^a^ if n is even 

Note that in both cases, we have 

with g smooth. Indeed, by using that we can write f'{p) = P^^^ f{p) with / smooth and positive, 
we have in both cases : 

g{a') = {a'y^~'{f{{ay))K (81) 

This is the natural generahzation of the change of unknown used in [2]. Then, thanks to the 
equation on a^, we get for {F ,u ) the system 

O^F'^ + ^0 . v^fc ^ 1 «o^(aO)V • u'' + V2g{a°) n^' • Va" + — f 1 + "°^/^f^ W • u' = V2g{a^)S':, 
V2 2 V g{a^) ) 

dtu^ + n° • Vm'^ + ^ V(a°5(a°), F'^) + u^ ■ Vn° = VS^. 
v2 

Note that the coefficient ° ^^o^ -^ is smooth even when aP vanishes since g is under the form ([8T]) . 
We have obtained a hnear symmetric hyperbohc system with a zero order term and a source term 
iS'^ depending only on (a-', (/j-') for < j < k under the form 

3=1 \ ^ J 

where A^(t, x) are smooth, real and symmetric and the matrix L is smooth. By the classical theory, 
there exists, on [0, T*), a smooth solution {F^,u^) in H°° x W^ of this system. Once u^ is built, 
we get a by solving the transport equation for a which is given by the first line of (j80p . Finally, 
we deduce the phase 99*^ by integrating in time the first line of ([78]) . We obtain 



/(i, x) = - f (2/'(|a0|2)(a°, a'=) + V/ • u^ - 5^) (r, x)dr. 



Finally, we choose in a similar way (a™^^, (/j™"*"^) that solve (j78|) with the initial condition 
(0™+^, c/?'""'"^) ,^^ = (oq, (^q) . Because of the assumption (j73|) . we find that they are also uniformly 

bounded in H^~'^ x H^ with respect to e. This concludes the proof of Lemma [2j D 

4 Nonlinear stability 

In this section, we give the proof of Theorem [2l We shall actually prove directly a more precise 
version which states the existence of a WKB expansion to any order. 
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Theorem 5 Consider \I'q = aQC^^o/'^ with Oq G H°°, c^g G H°° and that for some ?n, G N, there 
exists an expansion (I72p as in LemmalM We assume {A) and let {a^ , ip^) be the smooth approximate 
solution given by Lemma\^ which is smooth on [0,T*). Then, 

• if m = 0, there exists Eq > and T G (0, T*) such that for every e G (0, eo]? the solution of ([1]) 
with initial data ^q remains smooth on [0, T] and satisfies for every s G N, the estimate 



^'^ exp ( ip'^) 



< ae. 



L°^{[0,T],H=) 



• if m > 1, for every T G (0,T*), there exists So{T) > such that for every e G {0,eQ{T)], the 
solution of ([T]) with initial data ^q remains smooth on [0, T] and satisfies for every s G N, the 
estimate 



*^exp ( ip"^) - a^ 



< CsTE 



m+l 



L°°{[Q,T],H=) 

Note that Theorem [2] is actuahy the special case m = in Theorem [5l 

Proof of Theorem [5] 

Let s > d/2. We take (a^,(/3^) the approximate solutions given by Lemma [2] and look for the 
solution of ([T|) under the form ^^ = (a^ + w)e'^'^^ '^ . We get for w the equation (j2ip with F^ given 
by (j22p and the initial condition w/t=o = 0. For s > d/2, and every e > 0, this semilinear equation 
is locally well-posed in H^: we get very easily that there exists for some T"^ > a unique maximal 
solution w G C{[0,T^),H'^) of (I2ip (see [5j for example). We shall prove that T^ is bounded from 
below by some T > if ttt. = 0, and that T^ > T for every T G (0,T*) for e sufficiently small if 
m > 1. Let us define 

r" = sup{rG(0,T"), VtG[0,r], 2Nl{w{t)) < e2™+4}. 

Note that r^ > since w{{)) = and that by Sobolev embedding, we have, for t < t^, 

\\w{t)\\l^ < Kh-^NI{w{t)) < K2^2m+2 < ^2^ 

for some K independent of e. 

We will apply Theorem [3] with F^ given by ()22p . To estimate F^, we use the following lemma: 

Lemma 3 Let R > 0, s > d/2 and w such that \\w\\^^ < R, and F^ given by ([22I1 . Then, for a 
constant C depending only on ||a^(t)||,^s+2,oo o.i^d R, we have 



|^e||2 _|_ _L||Ijxi_F^II^ ^ /-<^2m+4 , /-,,2m ^^e 



\H' 



\H= 



-1 < Ce" 



Ce''^Nl{w) + C 



Nljw) fNIJw) 

A "T 



N!{w). 



We postpone the proof of Lemma [3] to the end of the section. We can first easily end the proof 
of Theorem [SJ Notice first that, by definition of ^", we have 



Ra = R'^ + ieAa' = Oh^^""^ ) + OH>^{e) = O^fc(e), 
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for every k, uniformly for < t < T, hence 

1 



Applying Theorem [3] and Lemma [3] with R = K, we infer that for < i < r^, 

jNl{w{t)) < C7e2'^+4 + Ce'^"'Nl{w{t)), 
which gives immediately, since wu^q = 0, that 

in the following cases: 

• for m = 0, 0<t<T with < T < T* sufficiently small independent of e, 

• for m> 1, T G (0, T*) is arbitrary, < t < T and e < SoiT) with eo(^) sufficiently small. 
As a consequence, r^ > T as desired and 

IPilL°°([0,T],/f«(R<*)) - '~'s,T£ 

It remains to prove Lemma [3l 

Proof of Lemma [3] 

We recall that F^ is given by 

F^ = R^ + Q{w) = R' + {a' + w) (f{\a' + u-p) - /(|a^P)) - 2{w, a')f'{\a'\'^)a'. 

As a first try, we could use the rough estimate 

Q^{w) = 0{\w\'^) as w^O, 

which would lead to 

mis + 4iiimQ^ii^._. < gihiii. < gA^fH^ 



which does not allow to conclude in the proof of Theorem [5] for m = and does not give a sharp 
result for the existence time if m = 1. To get the refined estimate of Lemma [3l the idea is then to 
use a Taylor expansion for Q^ w.r.t. w up to second order, and write 

Q^{w) = \w\^f'{\a'\'^)a' + 2f'{\a'\'^){w, a')w + 2a' f"{\a'f){w, a'f + G'{x, w), 

so that for fixed x, we have as if — > 0, 

G'{x,w) =0{\w\^). 

We turn now to estimate each term in F' . 
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Estimate for R" = ieK^ - R'^a". Thanks to dZlD, we have 

Moreover, since K^ is real- valued and since, from ()38p . Ima"^ = Ow'°^i^), we also have 



thanks to ((75]) . We have thus proven that 



e^ 



Estimate for G^{x,w). The estimate relies on Lemma [5] in the appendix. Indeed, it is clear from 
the Taylor formula that G^ may be written under the form 

(Reii;) hii{x,w{x)) + {Kew)(jia.w)hi2{x,w{x)) + {haw) h22{x,w{x)), 

where hn, hu, /i22 : IK"' x C ^ C are of class C°° and Vx G W^, hn{x,0) = hi2{x,0) = h22{x,0) = 
0. Moreover, hn, hi2 and /122 verify the hypothesis of Lemma [5] in the Appendix since o^ S 
L°°{[0,T],H'^). As a consequence, if ||'w||^oc ^ R^ 



II ii9 II 11*^ 

W\^, <C||u;||^,, 
which implies 

||G^(x,^x))||^. + l||lmG^(x,^x))||^._, < ^\\G'[xMx))fH^ < ^N^wf. 

The estimate for the quadratic terms in Q^{w) will rely crucially on the fact that o^ is real to 
first order and that {w,a^) is estimated in H'^~^ by Ng{w) and not just by e~'^N^{w). 

2 f'n^e\'2\^e 



Estimate for Ff = \w\ f {\a'^f)a'^ . We have 



\F!\\l^ < g^^|H^ 



and in view of (j38p . Im a"^ = 014/8,00(6), thus 



l||lmFf ||5^.„, < C\\\w\%,^, < ^Nliwf. 



Estimate for i<f = 2/ (|a^| ){w, a^)w. We begin with the rough estimate 
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Moreover, one has 

\\f'{\ar){w,a')\\;j,_,<CN!{w). (82) 

Indeed, let ^ € N*^ with \fi\ < s — 1. Then, 

d^^{f'i\a'\^)iw,a')) = Y. * d^[n\a^^)]{d''w,d^a'), 

where * is a coefficient depending only on a, (3 and A. Since |;u| < s — 1, the terms {d°'w, d^a'^) are 
bounded in L^ by S(?i;) 2 + e||Ti;||j:/s-2 as soon as \a\ < s — 2. The term in the sum with \a\ = s — 1 
(hence fi = a and /? = A = 0) is /'(|a'^p)(9'^zi;, a^) is bounded in L^ by N^{d^w). Hence, 
follows. 

As a consequence, by ([HT]) and Sobolev embedding, we obtain 



\\n\a^\'){w,anw\\^^_, < Cs\\w\\^^{\\n\a^\'){w,an\\Hs-. + \H\h^-^) ^ ^^IH- 
Consequently, 



|F|||^, + ^||lmF|||^,_, < ^NHwf. 



Estimate for F| = 2a^ f"{\a''\^){w,a''y. We find as for Ff 



\Fi\\l^ < ^.Nliwf, 



and once again in view of (j38 



We conclude the proof of Lemma [3] summing these estimates. D 

5 Geometric optics in a half-space 

In this section, we consider the Gross-Pitaevskii equation in a half-space in dimension d < 3 

GP(*^) = ie5t*^ + — A^^ - ^^(I^^P - 1) = 0, j;eM^ = M'^-^ X (0,+cx)). (83) 

We consider the Neumann boundary condition (J12p on the boundary and the condition (|13p at 
infinity, that is 

d^^ 9^^ / i r, i \ 

and exp \u°°\H - -u°° ■ x)"^^ ^ 1 |x| ^ +00 



dn /dR^ dz /z=o \ 2e e 

by using the notation x = {y, z) G M x (0, +cxd). 
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5.1 Construction of the WKB expansion 

In this section, we shall consider a smooth solution {a,u), with a real- valued, of 

dftt + u ■ Va H — a V • n = 

2 (84) 

dtu + u-Vu + V{a^) = 0, 
with the boundary condition Ud{t,y,0) = and the condition at infinity 

u{t,x) — > u°°, a{t,x) -^ 1 when |3;| -^ +cxd. 

Since we look for a real- valued, the resolution of this system is made in 1^14] (Theorem 2). Given 
s E N*, if the initial datum oq is positive and (oq — 1,uq — u°°) S H^, and under some compatibility 
conditions for (oq, ^^o) on the boundary dW^ of sufficiently high order on the initial data, there exists 
To G (0,+oo) andasolution (a,n) on [0,To] with {a-l,u-u°°) G C°{[0,To], H')nC'^{[0,To], H"''^), 
such that 

a{t, x)>a>0, Vt G [0, Tq], Vx G mJ . (85) 

for some a > 0. We also define the phase ip by 

(f{t,x)=^o{x)- I (^|Vv9|2 + |a|2-l)(r,x) dr. 

In view of the condition (|13p at infinity, (/? is not in H'^ but 99(t, .) — u^ ■ x — ||m°°|^ G H"^. 

The aim of this subsection, is to prove the existence of WKB expansion (which involves boundary 
layers since the solution of (I84p does not match the Neumann boundary condition (I12p ) up to 
arbitrary orders for (j83p . (J12p . (|13p starting from a smooth (a,ii) which verifies (j85p . 

We define the set of boundary layer profiles 5exp as 

5exp = {^(t,y,^)ei?°°(M+xM'^-ixM+), Vfc, a, /, 37 > 0, |5f 9^9^.4| < C^,,,, exp(-7Z)}. 

Lemma 4 Lei s G N 6e /ixerf. For every m £ N* , there exists a smooth approximate solution 
qfa,m _ Q£g*— Q^ [0,Tm] of (j83p . wii/i i/ie Neumann condition (|12p and i/ie condition (J13p ai 
infinity, such that 

GP{^^'"') = e'^R'e'^, (86) 

where R^ can he written under the form 

R' = -a^fe"*'"(t,x) + i?5^"(t,2/,-)) +i(ei?-*.-(t,x) + i?^''"(t,y,-)), (87) 



wzi/i i?<^ ' , Ra ' smooth and uniformly bounded in H'^ and Rd {t,y,Z), R^ (t,y,Z) G Sexp- 
Moreover, ^"'™ verifies (J12p . (J13p . a^ is real-valued and a^ , 99^ /lave smooth expansions under the 
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form 

m— 1 



fc=i 

m— 1 



^e ^ ^+Y,e''h\t,x) + <^\t,y,-))+e"'<^"'it,y,-). (89) 



k=l 
The boundary layer profiles A^{t,y, Z), ^^{t,y, Z) belong to Sexp o-nd are such that 
dzA^{t,y,0) = -5,a(t,y,0), dz'^Ht,y,0) = -a,^(i,y,0), 

dzA''{t,y,0) = -d,a''-\t,y,0), dz<^Ht,y,0) = -dz^''-\t,y,0) V2 < fc < m. (90) 

Proof : 

Since ^'^'"^ = a^ exp {i—), we want to solve approximately 

1 1 2 

Since, in this section, we are looking for a^ real-valued, we can split the system ()9ip into 



9ta^ + Vip^ • Va^ + - a^A93^ = 



for t > 0, X e Mi. (92) 



e 



2' '^ ' ' ' 2 a' 

Note that in this section, the division by a^ in the right-hand side of the second equation of (j92p is 
not a problem since a^ = a verifies (|85p and hence does not vanish. 

We thus plug the expansions ([88]) . ([89]) in ([92]) and we cancel the powers of e. To separate 
interior and boundary layer terms, we use the general theory of |JJJ. In particular, we use that for 
every function /, we have the expansion 

/ (u{t, x) + V{t, y, z/e)) = f (u{t, x)^ + f {u{t, y, 0) + V{t, y, z/e)^ - f (u{t, y, 0)) + eTZ, 

where TZ £ Sexp- This yields that the boundary layer part of f(u{t,x) + V{t,y,z/e)^ is given by 
f[u(t, y, 0) + V{t, y, z/e)) — f{u{t, y, 0)) . In the following, we use the notation Wb = W{t, y, 0) for 
every W{t,x). At first, the £~^ term in the equation only gives 

abdzz^^ = 

and hence we have <I>^ = 0, since a^ > a > and <I>^ G Sexp- Note that this is coherent with the 
fact that Ud{t,y,0) = dz^b = so that we do not need a boundary layer to correct the boundary 
condition. The e^ term gives 

dtf + ^\Vip\^ + a^-l = 

for t > 0, X e Rl (93) 

dta + Vip ■ Va + -a Aip = 
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for the interior part (which is the expected equation) and for the boundary layer terms, for (t, y) S 

TU>+ -^ 1U><i-l 



abdzz^^ = -{dzip)bdzA'=0 for Z > 0, 
since {dz(p)b = Ud{t,y,0) = 0. Consequently, we also find ^^ = 0. Next, the order e gives 



(94) 



dta^ + Vlp ■ Va^ + Vip^ ■ Va + -{aAip^ + a^Aip) = 

, dt^p^ + 2a a^ + Vv? • Vip^ = 
in the interior and for the boundary layer terms 
r 1 
< 
, abdzz^^ 



for t > 0, X £ 



\dzzA^ = A' (^dtif + ^\Vip\^ + a^ - l)^ + 2alA^ = 2alA^ 



for Z > 0, (95) 



G' 



where G^ G Sexp depends only on (0,^^,0"*^) and (99, c^^). Consequently, the boundary layer A^ is 
given by 



A^ 



2ab 



b_^~2atZ 



in order to match (j90h . Finally, the e , k > 2 term gives 



dta!" + V(f ■ Va!" + Va • Vip^ + -A^'' + —Ay? 



for t > 0, X G 



(96) 



s: 



and 



dzzA^ = AalA^ + F^ 



dzz<^' = G' 



for Z > 0, 



(97) 



where S^ and S^ depend only on {a,(p), (a-', (/j-')i<j<fc_i, F^ G Sexp depend only on {a,ip), 
{a^,ip^,A^,^^)i<j<k_i and {a^ , ip^ , ^'') , and G^ G 5exp depend on {a,ip), {a^ ,(p^ ,A^ ,^^)i<j<k-i. 
Therefore, if we want to solve by induction these equations, one has to determine first $ , then 
{a^,ip^) and finally A''. 

To solve the cascade of equations by induction, we first find {a^, (p^). As before, we notice that 
(a^,n^ = V(/?^) solves a symmetrizable hyperbolic system (there is no problem with the vacuum 
since we are in the same situation as in \9\). Since the condition at infinity is already absorbed by 
(a, (/?), one can look for {a^,u^) in H^. Moreover, we solve the system in M^ with the boundary 
condition u^{t,y,0) = which is needed in order to match (j90p since we have already found that 
<5^ = 0. The existence of a smooth solution for this linear system with the boundary condition 
u^{t,y,0) = which is maximal dissipative can be obtained by the classical theory [T^. Then, one 
finds 99^ by the formula 



ip^{t,x) 



(2a a + u ■ u )(t,x) dr. 
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Furthermore, since F'^ S Sexp and a^ > a > 0, the first equation in ([97|) (with k = 2) has a unique 
solution ^^ G Sexp- We have therefore found (a^, A^, t/?^, $\ A^, <l>^). 

We now proceed by induction. Assume that, for some m > 2, we have determined {a^ , if>')i<j<m-i 
and {A^ ,^^)i<j<rn- Then, we wish to solve (f96]l and (i97]) with k = m + I. Since G'""'"^ is already 



determined and G'"+^ € 5exp, the differential equation dzz^"^^^ = G™+^ has a unique solution in 



,m+i,_. ^^_ /+°°G-+i(t,y,C) 



'5exp and 

Jz ab{t,y) 

This determines the boundary condition for it™+^ = V(^™~*"^. Indeed, to match (|90p we shall need 
to impose 

n-+i(t,y,0) = (5,</.'"+i)(t,y,0) = -(az<I>™+i)(t,y,0) = / ^4^ '^C, (98) 

Jo afe(i,y) 

which is non-zero in general. We then solve (j97p in the following way: (a™~^^,ii™''*'^ = Vc^™^^) 
still solves a linear symmetrizable hyperbolic system, with source terms S™^^ and S^'^^ already 
known, with the maximal dissipative boundary condition ()98p . It has then a smooth solution by 
the above mentionned theory. Then, we recover cp"^^^ as usual by 



ip'^+\t,x)= / (5™+i-2aa"+i-'u-u™+i)(r,a;) dr. 



Finally, the first equation in (197p (with k = m + 1) is a linear ODE for A^^ , with source term 

^exp 



pm+1 g 5ea;p now determined, for which we can write down the unique explicit exponentially 



decreasing solution satisfying dzA''{t,y,0) = —dza^{t,y,0). 

Consequently, we have constructed an approximate solution of (j92p such that 

dta' + Vip' ■ Va^ + \a'A^' = e'"(i?f*''"(t,x) + e-ii?^''"(t,y, z/e)) 



dty,' + - iVv.^p + {a'f - 1 = ^^(t,x) +e-(i?-*''"(t,x) +i?i;™(t,y,z/e)) 



where RT '"*(*, x), i?™ '™(t,x) are smooth bounded functions and i?a'™, R(/^^ G ^e^^p. We can thus 
write the error i?^ in the GP equation as 

R'{t, x)=e"'(^- a'iR^''"" + R^f) + i(e<"*''" + i?^''")) . 

This ends the proof of Lemma [H D 

5.2 Validity of the WKB expansion 

We shall now prove the stability of the WKB expansion built in Lemma HI 

Theorem 6 Let "^"-'"^ = a^e* e a WKB expansion defined on [0,Tm] given by Lemma\^ Then for 
d < 3, and m> A there exists a unique smooth solution ^^ also defined on [0, T^] of (j83p . ()12p . 
(fT3l) such that '^%^q = ^/4=o' Moreover, we have the estimate 

e\\^'e^ - a'W^, + e^^l^^e"^^ - a'\\^, < C7„e™-5, Vt G [0,r^] 
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and in particular 

W^'e-'^ - {a + eA')\\^,^^ < C^max{e,e'"-i}. (99) 

Remark 3 For simplicity, we have restricted ourselves to dimension d < 3. Note however that it 
is possible to get H^ estimates for every s. By contrast with Theorem [2l we emphasize that the 
initial condition in Theorem [6] is exactly the WKB approximate solution ^i"-'"^. In particular, this 
initial datum has to verify some compatibility condition on the boundary. 

Proof. 

As in the proof of Theorem [SJ we set 

and we study the equation for w i.e. (fT9]l . Note that we are now seeking for w which tends to 
zero at infinity since the boundary condition at infinity is already absorbed in the WKB expansion. 
Again the first step is to get estimates for the linear equation (I2ip in O with the Neumann boundary 
condition 

a,w;(t,y,0)=0. (100) 

As we can check in the proof of Lemma [H in all the integration by parts that are performed, the 
boundary terms vanish due to the Neumann boundary condition or the fact that n^(t,y,0) = 0, 
and hence the proof of the L^ stability will be almost the same as the one in the whole space. 
Nevertheless, we have to pay attantion to the presence of boundary layer terms in the coefficients. 
At first, we note that since <I>^ = and <I>^ = in the WKB expansion, we still have that M (which 
is defined in Lemma [I]) is independent of e. Indeed, for the worse term which is V(V • u^), we have 

V(V • u') = dzzz-^^ + VAv? + 0{e). 

Next, keeping the definitions of Ra and R^ given in (fT7|) . (fT8|) and by construction of the WKB 
expansion, we have 

WRah^ < Ce^. (101) 

Nevertheless, again by construction of the WKB expansion, we only have 

_ e^/\a^ 

and due to the presence of boundary layers in a^ , we can split R^p into 

R^ = e^R^'\t, y, z) + ei?^(t, y, -^ (102) 

where i?™* is smooth and bounded whereas i?^ G S^xp and we see that e ||i?t,||/^oo = 0(e), 
e \\I1x^p\l°° = 0{1), hence (|23l) would be useless. Moreover, the fact that i?^ belongs to Sexp does 
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not seem to improve the estimates. The way to overcome this difficulty seems to incorporate this 
new singular term into the functional. Let us define the operator 

S^w = -—Aw + 2{w, a'')a^ + eR^^w, 
our weighted norm in this section will be 
Nliw) = I ({S%w,w)+Ke^\w\'^\dx = ]- I (e^\Vw\^ + A{w,a^f + 2eB'^\w\^ + 2Ke^\wf\dx. 

Note that R^p has no sign, nevertheless, N'^{w) can be bounded from below by a weighted H^ norm 
if K is chosen sufficiently large. Indeed, since i?^ belongs to Sexp we can write 



2e 



R^ \w\ dx 



<Ce I e ^^|w;pdx 



and then use the one-dimensional Sobolev inequality 

\w{t,y,z)\^ <C( [ \w{t,y,z)\''dzY ( [ |5,u;(t, y, z)p dz 
to get 

/ IX£. I 1 9 / 2L£. 9 

e e ^ \w\ < Ce\\w\\i^2 \\Vw\\i2 e ^ dz < Ce \w\12 \Vw\12. 
Jn Jn 



In particular, we have proven that 



2e 



R^ \w\ dx 



< Ce^ \w\l2 \\/w\i2. 



This yields thanks to the Young inequality 



2£ 



Rt\w\'^dx 



< -e'^WVwfri +C£'^\\wfr2 



(103) 



(104) 



(105) 



where C is independent of e. Consequently, if K is chosen such that 2K > C, we get 

Nliw) > Co(e^lklll^i + f {w,a'fdxy Co > 0. 

Note that in this section, we have 

a^ = a + 0{e) 

with a > a, this finally yields that N^{w) is equivalent to the weighted norm 

Nl{w) ~ e^ll^ll^i + ||Re u;||^2. (106) 

The first step in the proof of Theorem [6] is to prove the equivalent of Lemma [TJ We shall prove 
the estimate 

d 



dt 



iV=(u;(t)) < CN^{w{t)) 



(107) 



+Il^'lli2+ / -{w,a')iia',F' 
Jn £ 



{ieAw, F^ 



{iF',Kw) 



where C independent of e. 
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Proof of (fTUTD 

The proof follows the same lines as the proof of Lemma[TJ At first, since S^ is self adjoint, we have 

— / {Slw,w)dx= / (2{S^^w,dtw) + 4{w,a^){w,dta^) + 2edtRl\w\'^) dx. 
dt Jn Jn^ ' 

Since dtJx^p G Sex-p-, we can still use ()103p to get 

2e [ dtBi\w\'^ <CNl{w) 

Next, as in the proof of Lemma [U we use (|21|) to express dfW as 

i , 1 , e'^R^^ iF^ 

dtw = — Slw — (n^ • Vw + -wV ■ u^) — i ^—w 

e ^ 2 ' e e 

to get 

/■ /■ / 1 ie'^R^"'* F^ \ 

2/ (dtw,S'w)dx = 2 { - {u"" ■Vw + -wV -u") ^w - i — ,Slw)dx. (108) 

The only term in the right-hand side of (llOSp which is not present in ()25p . is 

J = -2 / ('u^ • Vu; + -u; V • u^ eR^^w 

Indeed, we have the cancellation 



[{iR^^'w,Rlw)=0 
Jn 



since i?"* and R'' are real. To estimate I, we note that we have a bound on the second term 
by using again (jl03p . It remains to estimate the first term. Integrating by parts and using that 
'u^(t,y,0) = 0, we get 



-2 [ {u' ■ Vw,£RIw) =e I V ■ {rIu') 
Jn Jn 



u;|2= / V ■ u' eRl\w\'^ + / eu' ■VR^^\w\'^. 
Jn Jn 



Again, the first term can be bounded thanks to ()103p . For the second one, we first notice that since 
ti^(t, y, 0) = and i?^ G Sexp, we have 

elu'-VR^^l <Ce(\VyRl\ + \zd,Rl\^ < Cee"?. 
This finally yields 

thanks to a new use of (jlOSp . 

The end of the proof of (fTUTD is then exactly the same as the proof of Lemma 1, since all 
the integration by parts do not create boundary terms either because of the Neumann boundary 
condition or because li^ vanishes on the boundary. 
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Higher order estimates 



The estimates of higher order derivatives are more involved than in the whole space. There are 
two main reasons. The first one is that there is a new singular term eR w which creates bad terms 
when we take the derivatives of the equation. The second reason is that to recover estimates on the 
normal derivatives, we need to use the equation which gives in particular that e^i9^ behaves like edt 
and eV. This anisotropy in the weights does not seem to allow to construct high order functionals 
like Ngi^w) which allows to get H'^ estimates without additional loss of e. Let us use the notation 



A = {Ao,--- Ad) = {dt,Vy,p{z)d,Y 

where the weight p{z) is given by p{z) = z/{l + z). Note that we can apply A to the equation 
since Aw still satisfies the Neumann boundary condition. The use of A is classical in hyperbolic 
characteristic initial boundary value problems (see [TTj for example) The weighted norm that we 
shall estimate is 

Yl{w) = Nl{w) + Nl{eAw). 

In dimension d < 3, this is sufficient to get the nonlinear stability. We shall see in the proof why 
the use of A^^ is necessary 
We shall prove that 

jYliw) < C (y^{w) + X%F') + X'ieAF')) (109) 

for some C > independent of e where we have set 

II PII2 llT-m Z?l|2 

^ [^) - W^ \\m + —^ + — ^2 • 

Proof of mrm 

As a preliminary, we shall rewrite pU7|) in a more convenient form. We can use that a^ = a-\-0{e) 
with a real, perform an integration by parts and use (jlOSp to get from (jl07p that 

jN'{w{t)) < CN'{w{t))+X'{F') (110) 

where 

II Z?l|2 llT-m Z7'£||2 

X [^ ) - jJ' 11^1 H — + ^2 • 

To prove (|109p . we start with the estimate of N^{edtw). When we apply edt to (f2T]) . we find 



iedt + Cjedtw = R^ edtw + edtF^ + C (111) 

where the commutator C can be splitted into 

C = Ci + C2 + C3 (112) 
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with 



Ci = edtR^w, 

C2 = 2e({dta^,w)a'' + {a'',w)dta' 



C3 = -ie'^ldtv!' ■ Vw + -dt{V ■ u'')w 

Consequently, we can apply (jllOp to (jllip with the new source term edtF^ + C to get 

j^Nl{edtw{t)) < CNl{edtw{t))+X'{edtF') + X'{C). (113) 

Thus it remains to estimate X^(C). Let us begin with X^(Ci). Thanks to the expansion (I102p . we 
easily get 



< Nliw). (114) 



Note that we could have a better estimate by using that R^ G Sexp and (|103p . Next, we turn to 
the estimate of X^(C2). By using that a^ = a + 0{e) with a real, we find 

X'iC2) < Nliw) + e||Re w\\l2 + e^llVt/;!!^ < N'iw). (115) 

Note that the above estimate was sharp. This is for the estimate of this commutator C2 that we 
had to chose the weight e in front of the time derivative. Finally, we estimate ^^(Cs). For the 
estimate of C3, we use that dtu^ vanishes on the boundary which implies that 

\dtu'd\<Piz). 

Thanks to this remark, we find 

X^Cs) < Nliw) + e^VAwWl, < Y^w). (116) 

Note that this is for the control of this commutator that we are obliged to add the vector field 
p{z)dz in the definition of the functional space. Consequently, the combination of (jllSp . (J114p . 
(fTT5]) and (fTTHI) gives 

^^Nl{edMt)) < Yl{w{t))+X^edtF^). (117) 

The estimate of eVyW follows exactly the same lines, and we also find 

±Nl{eVywit)) < Y^{wit))+X^{eVyF^). (118) 

The estimate of eAj^t;; = ep{z)dzW requires some additional work since the vector field A^ does not 
commute with the Laplacian. By applying eA^ to ([2T]) . we get 



iedt + C ) eAdW = R^eAdW + eA^F^ + C + C4 (119) 
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where C is defined as in ()112p above with dt replaced by A^^ and C4 is given by 



C4 = - — [Ad, A] = -—(2dzpd^zW + dzzpdzw). 



Next, we can apply (|107p to get 



d 



-NX {eKdw{t)) < Nl {eAdw{t)) + X%eAdF') + X%C) + \\Ci\\l, 
+- / (eAdW,a")(ia^C4)- / {iC4,BieAdw) 

Since one can easily check that X^{C) still satisfies the bounds (I114J1 . (|115p . (J116p . we obtain 
d 



dt 



Nl{£Adw{t)) < YX{w)+X%eAdF') + \\C4J,^ 



+- {eAdW,a'){ia',C4)- / {iC^, E^ eAaw) . 
£ Jn Jn 



Next, we note that 



and that 

4 

e 



{eAdW,a^){ia^,C4] 



<1 



£ Jn 



\\C4m ^ e'll«^llH3 



e|a,w;| |p(z)C4| < -2 Are 



1 



In a similar way, we also get 



{iCi,R'^eAdW 



Consequently, we have proven that 
d 



e Nl{w)2 ( \\pdzzw\\r2 + \\dzW 
£\\dzw\\L2\\pC4L'^ <YX{w). 

6||„,,||2 



dt 



NlisAMt)) < Yl{w)+X'{£AdF') + e^wfH,. 



(120) 



To conclude, it remains to estimate e^||w|||^3. As usual, this is done thanks to the equation ([21 
and the standard regularity result for elliptic equations. We rewrite (PT|) as the equation 



e'^Aw = G', dzwit,y,0) = 



(121) 



where the source term enjoys the estimates 



\m\l2<eMw\\l, + \\w\\l, + \\F%,, 
iVG^Ili^ < e'\\VAw\\l2 + Ml, + iVF^Ili^. 



Consequently, we get from (J12ip by standard elliptic regularity that 



(122) 
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By replacing this last estimate in (jl20p . we finally obtain 

j^Nl{eAMt)) < Yl(w) + X^eA^F^) + ||F^||^,. (123) 

To conclude, it suffices to sum the estimates (fTTOl) . (fTTTl) . (fTTSjl and (fT23]) to get (fT09]) . 

The estimate p09p is sufficient to prove the nonlinear stability stated in Theorem [6] for d <3. 
Nevertheless, it is possible to prove by induction that for every s, 

1(5] NliisArw)) < Y, {x^{ieArF^)+Nl{{eArw)). 

m<s m<s 

Nonlinear stability 

Thanks to ()109p and Gronwall inequality, we get for < T < Tm 

supyf (u;) < yf (0) + Te'^ sup (x^{F^) + X^{eAF^) 

[0,T] [0,T] ^ 

for some 7 > independent of e. Consequently, we can combine this last estimate with (jl22p to 
get 

sup Z%{w) < y^(0) + (1 + T)e^^ sup (x'{F') + X'ieAF')) (124) 

[0,T] [0,T] ^ ^ 

with 

Z^(u;)=y|(T/;)+e6||y;||^3. 

Thanks to this a priori estimate, one can easily prove by standard fixed point argument the existence 
of a unique solution of (J19p on some interval of time [0,T^] C [0, T^] such that Z^{w) remains 
finite. 

By using that if/t=o = and the equation to compute the time derivative, we find 

Yt{w{0)) < Coe^"". 
Moreover, using that F^ = e'^R^ + Q^j we have thanks to (|87p that 

(l + r™)e^^'" sup (X%{R') + XliAR')) < C(T„)e2™^i 

where Tm is the existence time of the approximate solution given by Lemma HI We can thus 
fix i? > Co + C{Tm) and define r"^ the maximal time for which the solution of (J19p satisfies 
Z^{w{t)) < i?e^™~^. As in the proof of Theorem[Sl we shall prove that for e sufficiently small, we 
have T^ = Tm- Thanks to (jl24p . we have for every T, T < t^, 

sup Z'^iw) < iJe^™-! + (1 + T)e^T ^^^ (x^Q^) + X'ieAQ')) . (125) 

[0,T] [0,T] ^ ^ 

Here, the expression of Q^{w) is given by 

Q%w) = a^'l'wl'^ + 2{a'',w)w + w\w\'^. 
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To conclude, we need to bound the right hand side of (jl25p . To estimate the nonhnear term, we 
use that for d < 3, we have 

Ml^<Wmw\\m, (126) 

which gives 

Hli^ < ^ < ^^'"-^ VtG[0,r^). (127) 

We shall take m such that 2m > 5 in order to get Iu^Iloo < 1 for t £ [0, r^). This implies 

mi'rn < (Iklli- + MU) M'm < ^^- 
Next, since H^{M.'^) C L'^ for d < 3, we also have 



£2 ~ £2 V- ' N-N7.-; ~ ^6 • 

Consequently, we have already proven that 

X^Qn < ^^- (128) 

Next, we evaluate X^{eAQ^). At first, we write 

^'llAgil^i < e'llA^^ll^i (Ikllioc + \\w\\U) + e2||At/;||i4 ||Vt/;||i4 (l + \\w\\l^) 
and by using for d < 3, the Sobolev embedding H^ C L^ and the Gagliardo-Nirenberg inequality 

iiv/iii4<ii/iii.iivviii, 

we get for < t < r"^: 

^2||A^£||2 ^ Z+{^) , ^2||v7„,,||2 ||„,,||l ||V27,,P < ^+0^ 

^ Im^-V II //I ;i -4 1^*^ II vu)||jyi||U)||^i II V wil^a ;::; -Q 
Finally, by similar arguments, we also have 

||gAQ^||2, ^ 1 1 ^||. ||2 II II 2 ^ Ki^f 

~2 ~ Im^^IIl4 Il^i'|li4 r:; P^wiij^i \\w\\j^i ^ -g 
We have thus proven that 



X%eAQn < ^i^. (129) 

Consequently, by using (fT25]) . (fT28]) . (fT29]) . we get 

sup Z^iw) < fe2m-i + C{R)il + T)e^^ sup ^±^. 

[0,T] [0,T] e 

By choosing m > 4:, this allows to get that t^ = Tm for e sufficiently small and that 

sup Z%{w) < Ce^™-!. 

lO,Tm] 

Finally, the estimate (j99p follows by Sobolev embedding. This ends the proof of Theorem [6l 
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A A Lemma about composition in Sobolev spaces 

During the proof of Lemma[3l we have used a result about composition in Sobolev spaces. This 
result is very standard when h does not depend on x (see, for instance, |18j). 

Lemma 5 Let R > 0, s £ N and h = h{x,w) £ C'*+^(M'^ x M^jM), satisfying h{x,0) = for all 
X G M . Assume moreover 

A = snp{\\d^d^h\\^^^^,^^^y aGN^ /3 G N^, |a| + |/?| < s + l} < +cx,. 

Then, there exists C , depending only on A, s and R, such that, for any w G H''^{M.'^) satisfying 
\w\^cx,nid\ < R, we have h{x,w{x)) G H'^iW^) and 



\h{x,w{x))\^, < C\w\^,. 

Proof. The proof is by induction on s G N and relies on the Gagliardo-Nirenberg inequality. If 
s = 0, it suffices to notice that since /i(x, 0) = 0, then for w G Br., 

\h{x,w)\ < A\w\. 

Assume then the result for s — 1 G N. Let ^u G N"^ with |/i| = s. One has easily 

d^'{h{x,w{x))) =Y,*{dxdi^^h){x,w{x)){d^wiY {d^W2)\ 

where a G N*^, a < /x, /3, 7 G N^, p, g G N* depend on (i and 7, \a\ + p\(5\ + q\'~^\ = s, and * is 
a coefficient depending only on //, a, [3 and 7. Furthermore, since w G H^ n L°°, the Gagliardo- 
Nirenberg inequality yields, for 1 < A; < s, 



k 

,2s <Cfc,,||Ti;||;^,||u;||^: 






As a consequence, by interpolation, if u; G H'^ n L°° and ||ti'|| r^o ^ R^ then for 7 G N"', I7I < s, and 

2 

II o'Y II ^ /^ II Up 

||C'w||^p < Cs,p,r\\w\\'^,. 

Therefore, in view of |a| -|-p|/3| + q\'^\ = s, by Holder inequality, we can estimate the terms in 
d^{h{x,w{x))) for which a^ ^ (thus \a\ < s) as 

\\{d^dl,+^h){xMx)){d^w^Y{d'^W2Y\\^, < ^||aVf _|.| ||9^«^2i' .-iH < Cs,p,rAMhs- 

L IPT L TiT 

For the term for which a = /i, we note that since h{x,Q) = for x G W^, then (5"/i)(2;,0) = for 
any x G M'', so that \iw £ BrC M^, 

\{d'^h){x,w)\<A\w\, 
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which iniphes 

\\{d"h){x,w{x))\\j^2 < ^1^11^2 <^||^i'||^s 

Combining these two estimates gives 

\\d''{h{x,w{x)))\\^^ < Cs,p,rA\\w\\^, 
and the proof of the Lemma is complete. D 
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